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1 12.  SPONSORING  MILITARY  ACTIVITY 


The  techniques  of  Sea  Photo  Analysis  have  been  used  to  obtain  surface  truth  descrip¬ 
tors  of  the  open  ccean  from  an  aircraft.  Aerial  photographs  were  taken  such  that  the 
density  fluctuations  of  the  negatives  were  related  to  surface  slope.  The  slope  spectra, 
obtained  by  Fourier  analysis  of  the  negatives,  were  fitted  to-ari  equivalent  form  of  the 
Pierson-Moskowitz  wave  height  spectrum.  The  photographs,  the  wave  height  spectrum,  and 
the  rms  slope  equations  determined  by  Cox  and  Muuk  then  yielded  the  averaged  heading 
of  the  waves,  the  equivalent  wind  speed,  the  rms  wave  height,  and  the  rms  slope  for  dean 
and  slick  surfaces.  The  optical  techniques  of  Sea  Photo  Analysis  used  in  this  program  and 
digital  analysis,  used  when  whitecaps  were  present  in  the  photographs,  are  discussed.  The 
empirical  equations  used  to  determine  rms  wave  height  and  slope  and  auxiliary  rms  si  >pe 
and  wave  height  spectrum  equations  are  given.  Surface  descriptors  were  inferred  for  seas 
driven  by  winds  from  3  to  21  knots  occurring  on  16  separate  days. 
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ABSTRACT 


The  techniques  of  Sea  Photo  Analysis  have  been  used  to  obtain 
surface  truth  descriptors  of  the  open  ocean  from  an  aircraft.  Aerial 
photographs  were  taken  such  that  the  density  fluctuations  of  the 
negatives  were  related  to  surface  slope.  The  slope  spectra,  obtained 
by  Fourier  analysis  of  the  negatives,  were  fitted  to  an  equivalent  form 
of  the  Pierson-Moskowitz  wave  height  spectrum.  The  photographs, 
the  wave  height  spectrum,  and  the  rms  slope  equations  determined 
by  Cox  and  Munk  then  yielded  the  averaged  heading  of  the  waves, 
the  equivalent  wind  speed,  the  rms  wave  height,  and  therms  slope  for 
clean  and  slick  surfaces.  The  optical  techniques  of  Sea  Photo  Analysis 
used  in  this  program  and  digital  analysis,  used  when  white  caps  were 
present  in  the  photographs,  are  discussed.  The  empirical  equations 
used  to  determine  rms  wave  height  and  slope  and  auxiliary  rms  slope 
and  wave  height  spectrum  equations  are  given.  Surface  descriptors 
were  inferred  for  seas  driven  by  winds  from  3  to  21  knots  occurring 
on  16  separate  days. 
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DETERMINATION  OF  OCEAN  SURFACE  DESCRIPTORS  USING 
SEA  PHOTO  ANALYSIS  TECHNIQUES 


INTRODUCTION 

The  Naval  Research  Laboratory  has  been  providing  interpretation  and  quantitative 
analysis  of  aerial  photography  of  the  ocean  surface,  and  related  consultive  assistance,  to 
the  Federal  Aviation  Administration,  Department  of  Transportation.  The  analysis  is  being 
used  by  FAA  in  its  study  of  multipath  effects  in  communications  between  satellites  and 
aircraft  flying  over  the  ocean.  The  techniques  of  Sea  Photo  Analysis  [1] ,  developed  at 
NRL,  have  been  used  to  obtain  the  ground  truth  descriptors  of  the  ocean  necessary  in 
attempting  to  correlate  multipath  effects  with  the  roughness  of  the  scattering  surface. 

The  ocean  surface  roughness  may  be  described  in  terms  of  two  components— sea  and 
swell.  Sea  is  the  roughness  generated  by  the  local  wind.  light  winds  ripple  the  ocean 
surface  with  short-wavelength,  small-amplitude  wavelets  which  decay  rapidly.  If  the  wind 
persists,  the  wavelets  are  constantly  regenerated  and  longer  wavelengths  appear.  As  the 
wind  increases  and  persists  over  a  sufficiently  large  area,  the  proportion  of  the  longer 
wavelength  waves  increases  and  the  ocean  becomes  rougher.  Therefore  all  wavelengths, 
with  small  probabilities  of  existence  for  the  extremely  long  and  short  waves,  constitute 
the  roughness  referred  to  as  sea.  Swell  is  the  remnant  of  a  sea  generated  in  the  area  by  a 
wind  field  that  existed  at  some  previous  time  or  that  propagated  from  a  sea  generated 
elsewhere,  possibly  thousands  of  miles  away.  As  time  passes  since  the  wind  field  has  died 
out  or  moved  on,  or  as  the  swell  system  propagates  out  of  a  generating  region,  the  short 
wavelengths  disappear  first,  with  the  successively  longer  wavelengths  decaying  at  slower 
rates.  Thus  swell  commonly  consists  only  of  the  long  wavelengths  generated  by  a  wind. 

The-Naval  Research  Laboratory  has  been  analyzing  aerial  photographs  to  determine 
the  parameters  describing  the  sea.  This  type  of  analysis  is  generally  not  applicable  to 
swell  but  is  limited  to  the  relatively  short  wavelengths  present  in  a  sea  by  the  clarity  of 
the  earth’s  atmosphere,  the  stability  of  airborne  platforms,  the  quality  of  presently  avail¬ 
able  optics,  and  the  physical  and  chemical  properties  of  photogiaphic  film.  The  FAA  has 
been  obtaining  supplementary  data  on  swell  through  hindcasts  by  the  Naval  Oceano¬ 
graphic  Office  (see  the  Appendix), 

In  Sea  Photo  Analysis  (SPA)  the  physical  processes  by  which  the  ocean  surface  is 
depicted  in  a  negative  are  identified  to  yield  the  relationship  between  the  surface  parame¬ 
ters  and  the  optical  density  of  the  negative.  In  the  type  of  analysis  used  here  the  density 
is  given  in  terms  of  the  surface  slope.  When  the  photographic  negative  is  optically  an¬ 
alyzed,  the  intensity  of  the  diffraction  pattern  is  proportional  to  the  slope  spectrum  of 
the  ocean.  The  valid  portion  of  the  spectrum,  which  is  restricted  by  the  size  of  the  area 
that  can  be  photographed  and  the  resolution  of  the  photographic  process,  is  fitted  to  an 
accepted  form  of  the  spectrum  encompassing  all  wavelengths.  The  total  slope  spectrum 
together  with  other  experimentally  determined  relationships,  yield  the  surface  descriptors. 
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The  next  three  sections  of  this  report  are  devoted-to  the  mathematical  basis  cf  this 
type  of  analysis.  The  relationship  between  the  density  of  the  negativ.  and  the  surface 
slopes  of  the  ocean  is  developed,  and  it  is  shown  that  optical- analysis  yields  the  slope 
spectrum.  A  slope  spectrum  corresponding  to  the  optical  spectrum  is  derived  from  the 
Pierson-Moskowitz  wave-height  spectrum,  and  the  empirical  relationships  used  to  deter¬ 
mine  the  rms  wave  height  and  rms  slope  are  given. 

The  theory  is  then  extended  to  computer  analysis  of  the  photography,  which  had  to 
be  used  in  eight  of  the  sixteen  sets  of  data  that  could  be  analyzed.  Next,  minor  distor¬ 
tions  in  the  optical  spectrum  and  ways  to  avoid  them  are  described.  Finally,  the  data  and 
results  of  the  data  analysis  are  discussed. 


RECORDING  SURFACE  SLOPES 

The  local  coordinate  system  used  to  develor-  the  theoretical  basis  of  the  analysis  is 
shown  in  Fig.  1.  In  this  particular  study  the  ocean  was  photographed  from  an  aircraft 


z 
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Fig.  3— The  local  cpordinate  system  used  in  Sea 
Photo  Analysis 


flying  in  the  minus-X  direction,  with  a  possible  Y  component  dependent  on  analysis 
criteria.  The  camera  was  positioned  such  that  its  optical  axis  was  pointing  aft  and  down¬ 
ward  at  an  angle  of  40°.  The  position  of  the  origin  of  the  coordinate  system  shown  in 
Fig.  1  is  variable  and  located  at  the.  central  point  of  the  area  being  analyzed.  In  most  of 
the  analyses  the  aircraft  would  be  flying  directly  in  the  minus-X  direction  with  the  unit 
vector  C,  in  tho  XZ  plane,  directed  toward  the  camera.  The  variable  in  the  location  of 
the  origin  would  then  be  the  angle  0,  with  6  being  equal  to  40°  in  the  center  of  the 
negative. 

If  a  single  sinusoid  were  propagating  on  the  ocean  surface  at  an  angle  ^  with  respect 
to  the  X  axis,  the  local  coordinate  system  would  appear  as  shown  in  Fig.  2.  At  the 
origin  the  steepest  gradient  on  the  water  surface  makes  an  angle  <j>  with  respect  to  the 
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Fig.  2— The  local  coordinate  system  for  a  single  sinusoid  propa¬ 
gating  on  the  ocean  surface  in  a  direction  'k  with  respect  to  the 
Xaxis 


XY  plane.  The  normal  to  the  water  surface  at  that  point  is  represented  by  the  unit 
vector  N,  making  an  angle  4>  with  respect  to  the  Z  axis.  In  Fig.  2  the  angles  <p  are  greatly 
exaggerated.  The  angles  generally  encountered  in  this  type  of  analysis  are  less  than  15°. 
The  maximum  rms  slope  determined  in  this  study  was  14°,  with  an  average  of  10°.  Al¬ 
though  steep  slopes  do  exist  on  the  ocean  surface,  the  majority  of  such  slopes  are 
associated  with  wavelengths  shorter  than  those  analyzed  here.  The  steep  slopes  that  do 
enter  into  the  analysis  occur  infrequently  and  degrade  the  results  only  slightly.  The 
rational  for  this  will  be  explained  in  the  section  on  Minor  Spectral  Distortions. 

The  light  that  is  specularly  reflected  toward  the  camera  from  a  facet  of  water  surface 
comes  from  a  direction  represented  by  the  unit  vector  £.  pie  three  vectors  C,  N,  and  S 
are  coplanar,  with  the  angle  between  C  and  N,  and  §  and  N,  being  0.  The  variation  cf  0, 
and  hence  of  the  intensity  reflection  coefficient  of  the  water  surface,  with  slope  angle  4> 
is  one  of  the  factors  enabling  the  recording  of  surface  slope  angles. 

Figure  3  shows  the  variation,  in  the  intensity  reflection  coefficient  with  0  for  hori¬ 
zontally  polarized  light.'  Thus  the  light  from  the  sky  that  is  reflected  toward  the  camera 
is  modulated  by  the  slope  angle  of  the  water  surface.  A  second  factor,  which  would  in 
theory  enable  one  to  record  surface  slope  angles,  is  the  usual,  under  clear  sky  or  uni¬ 
formly  overcast  conditions,  monotonic  variation  of  sky  brightness  with  tj.  'Die  sky  is 
usually  brighter  at  the  horizon  than  at  the  zenith.  Since  both  the  sky  brightness  and 
intensity  reflection  coefficient,  in  crease  for  positive  q>  (away  from  the  camera)  and  de¬ 
crease  for  negative  <p  (toward  the  camera),  the  contrast  between  different  slope  angles  is 
increased.  The  presence  of  broken  clouds  may  degrade  the  analysis  somewhat,  that 
degradation  increasing  with  the  brightness  fluctuations  and  the  spatial  scale  of  those 
fluctuations. 

The  light  directed  toward  the  camera  from  a  given  point  on  the  surface  and  from- 
the  intervening  atmosphere  may  be  expressed  by 


/  =  bLS( 77,  $)R(P)  +  b2W(C,  0)  +  h B  , 
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Fig,  3— The  variation,  in-the  intensity  reflection  coefficient  R{(i) 
with  P  for  horizontally  polarized  light 


where  I  is  the  intensity  of  the  light  striking  a. corresponding  point  on  the  film,  bi,  &?,, 
and  63  are  constants  that  need  not  be  determined  in  a  first-order  development  for  the 
open  ocean,  S(r?,  f)  is  the.  brightness  of  the  sky  in  a  direction  determined  by  17. and  ?, 
R(fi)  is  the  intensity  reflection  coefficient  for  the. slope  reflecting  S(q,  f)  toward  the 
camera,  W(C,  4>)  is  the  intensity  of  the  light  scattered  upward  through  the  surface  making 
an  angle  <p  with  the  horizontal  and  in  the  direction  5,  and  5  is  the  intensity  of  light 
scattered  by  the  intervening  atmosphere  into  the  camera*  The  third  term,  63S,  is  rela¬ 
tively  uniform  over  the  field  of  view  of  the  camera,  possibly  increasing  toward  the  far 
field.  It  is  however  not  a  function  of  <p;  Also,  the  conditions  under  which  one  may 
obtain  good  photography  are  such  that  the  third  term  is  much  smaller  than  the  first. 
Since  the  third  term  is  not  a  function  of  <p  and  its  magnitude  is  very  small,  its  effect  on 
the  development  that  follows  will  be  ignored.  The  second  term  is  a  weak  function  of  <p 
but  may  have  a  significant  magnitude.  Its  effect  is  greatest  in  turbid  waters,  where  the 
relative  intensity  of  the  light  scattered  from  beneath  the  water  surface  and  refracted  to¬ 
ward  the  camera  is  large.  Even  in  this  case  the  variation  in  W(C,  <p)  with  4>  is  small  com¬ 
pared  to  the  variation  of  the  first  term  with  <p.  For  the  nonturbid  open  ocean  the 
variation  of  W{C,  <j>)  with  <P  is  very  small  compared  to  the  variation  of  the  first  term 
with  <p. 

To  a  first-order  approximation  the  intensity  due  to  a  given  slope  angle  is 


/  =  /n  + 


I-C 

<1 


¥: 


t£ 
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where  the  subscript  zero  indicates  that  the  term  is  to  be  evaluated  at  the  analysis  angle  0 

with =  0.  A  first-order  approximation  is  valid  because  the  slope  angles  being  analyzed 
are  small. 


=b  s  —  %L  bRdS  W 

d<t>  blS  dp  a 0  bl^.d4>  a 

^  b\S  cos1^  ~+  2b  iR  cos  ^  —  +6o  , 

op  A  3jj  *  30  ’ 


(3) 

(4) 


assuming  that 


as 

ar 


0.  Hence 


/  -/0(1  +  U0  <t>  cos  S')  , 
where  the  sensitivity  function  U0  is 

(dR\  .  „  „  fdS\  02 


^0  = 


f)o 


(5) 


(6) 


fhe  third  term  in  Eq.  (6)  is  very  small  for  the  open  ocean  and  may  be  ignored.  The 
practical  value  for  U0  is  therefore 


Rq  \90/o-  S0  \dvl  o 
This  value  is  nominally  0.1  per  degree. 

If  the  sea  is  photographed  on  the  linear  portion  of  the  7  curve  ( D  vs  log  E  curve) 
of  the  film,  the  optical  density  of  the  negative  is  related  to  <p  by 


(7) 


D  ~D0  +ylog{l  +  U0  <p  cos'if)  . 


(8) 


The  above  development  has  been  carried  out  for  one  specific  point  on  the  negative. 
Equation  (8)  is  however  an  equally  good  approximation  fqr  surrounding  points.  It  may 
therefore  be  written 


D{x,  y)  =  D0(x,  y)+  7  log  [1  +  U0(x,  y)  <j>(x,  y)  cos  *(*,  y )] 
=  Po(x,  y)  +  7  log  [l  +  <p\x,  y)]  , 
where  x  and  y  are  coordinates  on  the  negative  and 

4>{x,  y)  -  U0(x,  y)  <p{x,  y)  cos  #(x,  y)  . 


(9) 

(10) 

(11) 


e 


ROGER  O.  PILON 


OPTICAL  AN  ALYSIS 

The  negative  that  is  obtained  by  photographing  the  sea  is  th&efore  a  two-dimensional 
record  of  the  random  surface  slopes  that  are  generated  by  the  wind  and  any  swell  that  may 
have  been  present.  The  most  direct  method  of  analyzing  such  a  two-dimensional  record 
is  optical  analysis. 


The  system  used  in  this  study  is  outlined  in  Fig.  4.  The  He-Ne  laser  emits  red  ' 
light  at  a  wavelength  of  6328  A.  The  lens  in  the  spatial  filter  focuses  the  beam  through  a 
pinhole.  The  region  around  the  pinhole  blocks  any  nonparallel  rays  emitted  by  the  laser 


NEGATIVE 


which  would  appear  in  the  transform  plane  as  noise.  The  transform  lens  then  focuses  the 
light  to  a  point  at  the  transform  plane.  When  a  negative  is  inserted  into  the  converging 
beam,  the  density  fluctuations  representing  the  random  angles  of  surface  slope  diffract 
thv.  light,  giving  rise  to  a  Fresnel  diffraction  pattern  in  a  central  region  about  the  point 
imaged  by  the  transform  lens.  The  intensity  pattern,  a  sample  of  which  is  shown  in 
Fig;  5,  is  given  analytically  by  [2] 


If{xf,  yf)  =  (XF)-2 


J  A(x,  y)e-/<2 *imx*fVyf)  dxdy 


(12) 


where  Xf  and  yf  are  coordinates  in  the  transform  plane,  x  and  y  are  coordinates  on  the 
negative,  X  is  the  wavelength  of  the  laser  light,  F  is  the  distance  from  the  negative  to  the 
transform  plane,  and  A'(x,  y)  is  the  amplitude  of  the  light  immediately  after  passing 
through  the  negative.  The  pattern,  as  shown  in  Fig.  5,  is  symmetric  about  the  origin, 
with  each  of  the  two  lobes  containing  the  same  information. 

The  integral  in  Eq.  (12)  is  the  Fourier  transform  J  of  A(x,  y).  Hence 


If(Xf,yf)  =  (hFtm{A(x,y)]  2. 


(13) 
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The  intensity  distribution  in  the  transform  plane  is  then  proportional  to  the  Fourier  trans¬ 
form  squared,  or  spectrum,  of  the  light  amplitude  transmitted  by  the  negative.  It  wili  be 
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Fig.  5— A  representative  diffraction  patter.--slope  spectrum  from  I 

one  negative.  The  dashed  line,  orientedVat .  an  angle  *^It>  with  , 

respect  to  the  xf  axis,  is  the  average  line  of  symmetry'  of  the  j 

intensity  distribution. 
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shown  that,  in  the  region  of  the  transform  plane  to  be  analyzed,  the  intenvxy  distribu¬ 
tion  is  proportional  to  the  slope  spectrum  of  the  sea.  j 

1 

i 

The  amplitude  of  The  light  in  a  laser  beam  is  a  gaussian  function, and  may  be  J 

written 

f  : 

atea-a*-**'**4*4, 

! 

where  ^40  is  the  amplitude  in  the  center  of  the  beam,  x  and  y  are  measured  perpendicular  j 

to,  and  from  the  axis  of,  the  beam,  and  a  is  the  distance  from  the  axis  at.  which  the  j 

amplitude  is  1/e  of  that  on  the  axis.  Immediately  upon  passing  through  the  negative  the 
beam  has  the  amplitude  \ 

t 

A(x,  y).=  A0e^x2+y2)/a2lQ-D/x'y)/2  (15)  j . 

=  A  oe"(^2+^)/a210-Do(x)y)/2[1  +y(jc>  y)]  -y/2  ,  (16)  | 

Dg(x,  y)  is  a  monotonic  function  which  increases  from  nea~  to  far  field.  The  spatial  re-  | 

lationships  at.  the  negative  are  such  that  the  term  outside  the  brackets  is  a  very  slightly  \ 

deformed  gaussian  beam.  Denoting  this  term  by  g '  and  in  keeping  with  a  first-order  | 

development,  Eq.  (16)  becomes  f 

$ 

>Y 

Mx,  y)  *  g '  -‘V&7 g  '<i>'{x,  y) .  (17) 


$ 
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Consequently  Eq.  (13)  may  be  written 

If{xf,  yf)  =  (XF)-2  j|G'|2  -  yRe[G’(G'  ©  /)*]  +  'Ay 2|G*  0  ^|2  j  ,  (18) 

where  if  ig)  =  G'  J  ($')  =  Re[  1  indicates  the  real  part  of  the  quantity  in  brackets, 
the  symbol  ©  denotes  a  convolution,  arid  the  symbol  *  denotes. the  complex  conjugate. 

The  first  term  in  the  braces  is  centered  at  the  origin  and  is  referred  to  as  the  dc 
term.  It  is  the  Fourier  transform  of'the  distaribution  of  light  that  impinges  on  the  nega¬ 
tive  after  being  slightly  deformed  by  a  function  of  Dq(x,  y).  If  the  gaussian  were  not 
slightly  deformed,  |G'|?  would  itself  be  gaussian  iri  shape  and,  physically,  the  image  of 
the  gaussian  distributic*’  of  light  passing  through  the  pinhole  of  the  spatial  filter.  The 
effect  of  Do(x,  y)  is  to  very  slightly  broaden  the  dc  term.  Nearly  all  of  the  light  energy 
transriiittec.  by  the  negative  is  focused, into  this  term,  giving  it  a  peak  intensity  lOOO  or 
more  times  greater  than  that  of  the,  second  and  third  terms  combined.  It  is,  however, 
essentially  zero  a  short  distance  from  the  origin.  Its  effect  is  inconsequential  in  the 
method  of  analysis  used  here. 

The  second  term  is  weighted  by  G  'and  therefore,  in  the  region  outside  the  dc  term, 
is  small  compared  to  the  third  term  in  the  squar*  brackets. 

The  third  term  gives  rise  to  thn' bow-tie-shaped  pattern  shown  in  Fig.  5.  The  effect 
of  the-convolution  is  such  that  each  point  in  the  intensity  distribution  is  riota  “point” 
intensity,  but  rather  a  “spot”  having  the  same  shape  as  the  dc  term.  If  the  distribution 
is  sampled  using  a  detecting  area  larger  than  each  convolved  spot,  Eq.  (18)  becomes 


yf)  =  (constant)  \4>  'I2  , 


(19) 


where  the  prime  on  /^indicates  the  distribution  outside  the  dc  term,  when  sampled  by 
an  appropriately  sized  detector. 

The  term  Uq(x,  y)  in  Eq.  (11)  is  very  similar  in  its  effects  on  the  diffraction  pattern 
as  Dq{x,  y)  in  Eq.  (16).  Uq(x,  yvis  a  montonically  increasing  function  from  near  to  far 
range.  Writing 


0 y  [U0(x,  y)  4>{x,  y)  cos  ❖(*,  y)]  (20) 

=  5  [U0(x,  y))  ©  IT  [0(x,  y)  cos  #(*,  y)]  ,  (21) 


it  is  seen  that  Uq(x,  y)  slightly  broadens  each  convolved  spot  in  the  diffraction  pattern. 
Equation  (19)  then  becomes 


yf)  =  (constant)!  J[0(x,  y)  cos  'F(.r,  y))  |2 


(22) 


i.-  jLi.i .  j .’,'4L 
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With  reference  to  Fig.  5,  the  intensity  distribution  along  arty  radial  line  at  angle 
arises  from  all  slope  angles  having  approximately  the  same  'i'.  Due  to  perspective  the 
angle  is  hot  equal  to  \P,  but  is  given  by 


tan  it  = 


tan  ^ 
sin  0 


(23) 


The  dashed ‘line,  at  an  angle  with  respect  to  the  Xf  axis,  is  the  average  line  of  sym¬ 
metry  of  the  intensity  distribution.  If  the  sea  is  photographed  such.that4he  X  direction 
is  somewhat  parallel  to  the  wind  direction  the.i  ipw  represents  ❖k,. 


In  the  type  of  analysis  used  in  this  study  the  diffraction  pattern  of  each  photograph 
was  read  along  the  line  defined  by  \j/w.  Different  points  on  the  negative,  however,  con¬ 
tribute  to  the  intensity  distribution  alortg  these  lines  with  slightly  different  values  of 
\P(x,  y).  The  effective  deviation  in  ^(x,  y)  is  small,  less  tiian  ±  2°,  and  cos  T'fx,  y)  may 
be  considered  to  be  a  constant.  Equation-  (22)  then  becomes 


If(r,  i>w)  =  (constant)!  J [0(x,  y)]|?  , 


(24) 


where  lUr,  \jjw)  is  the  intensity  along  the  line  defined  by  at  a  distance  r  from  the 

j/iC^vAvivudl  tv't/iic  navel  TTavcieiigi/iu'-n^;  uu  uiC  Du:aii 


... 


and  may  be  related  to  k(=2n /\w)  when  the  altitude  of  the  aircraft  and;  certain  camera  and 
optical  bench  parameters  are  known.  Since  the  surface  slope  angles  being  analyzed  are 
small,  the  angle  is  essentially  equal  to  the  surface  slope.  The  intensity  distribution  may 
then  be  written  in  a  form  more  familiar  to  the  oceanographer: 


I'fiK  i//w)  =. (constant)  <l>(0;  k,  0) ,  (25) 

where  <I»(0;  k,  0)  is  the  slope  spectrum  of  the  sea  surface,  given  as  a  function  of  k  for 
slopes  in  the  direction  of  the  wind  (^  -  'I'm  =  0).  An  example  of  the  distribution  If  is 
shown  in  Fig.  6. 

In  the  preceeding  development  it  has  been  assumed  that. the  area  of  the  detector 
sampling  the  intensity  distribution  is  larger  than  any  convolved  “spot.”  This  is  not  a 
limitation  in  the  type  of  analysis  used  here.  Any  given  negative  depicts  only  a  small 
sample  of  the  random  slopes  on  the  ocean  surface.  The  inten^y  distribution  is  therefore 
not  uniform,  but  a  pattern  ox  “spots”  representing  only  the  wavelengths  in  the  area  being 
analyzed.  The  detecting  area  must  be  large  enough  to  integrate  the  distribution  to  yield 
a  reasonably  smooth  curve,  such  as  shown  in  Fig.  6. 


THE  OCEAN  SPECTRUM 

The  ideal,  way  to  analyze  the  ocean  surface,  and  obtain  such  parameters  as  rms  wave 
height  and  rrns  slope,  would  be  to  obtain  one  or  more  negatives  depicting  an  area  suf¬ 
ficiently  large  so  that  the  longest  water  wavelength  present  would  be  shown  m  my  times, 
and  at  the  same  time  adequately  resolve  the  shortest  wavelengths  present.  Since  all  slopes 
present  would  then  be  represented  many  times,  a  more  accurate  spectrum  could  be 
determined  if  the- negative  showed  no  perspective.  Obtaining  such  a  negative  is  clearly 
not  possible.  One  could  possibly  construct  the  spectrum  from  a  number  of  negatives, 
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Fig.  6— A  representative  example  of  I  f(>,  tyw)  vs  r 


each  recording  a  different  wavelength  portion  of  the  ocean  surface,  but  such  an  under¬ 
taking  would  be  extremely  difficult  in  practice. 

In  the  type  of  analyses  used  here  it  is  necessary  to  record  only  one  band  of  wave¬ 
lengths  from  roughly  2  to  50  meters.  The  shape  of  the  optically  determined  spectrum  is 
then  compared  to  a  previously  arid  experimentally  determined  ocean  spectrum  in  the  same 
wavelength  band.  Since  only  the  shape  of  the  optical  spectrum  is  used,  it  is  not  necessary 
to  determine  the-constant  in  Eq.  (25).  This  method  effectively  extends  the  limited 
optical  spectrum  to  all  wavelengths,  thus  yielding  a  more  accurate  determination  of  the 
ocean  surface  descriptors. 

The  ocean  spectrum  to  which  the  optical  spectra  were  fitted  is  that  determined  by 
W.J.  Pierson,  Jr.,  and  L.  Moskowitz  [3].  The  wave-height  spectrum  as  a  function  of 
frequency  is  given  as 


,  co)  =  8.1  x  10~3 


-fe4/v4w4 


where  h  is  the  wave  height  (displacement  from  average  surface  level),  c o  is  the  temporal 
frequency  of  the  wave,  g  is  the  acceleration  due  to  gravity  (32  ft  s-2),  0  =  0.74,  and  V  is 
the  wind  speed  measured  64  feet  above  the  mean  ocean  surface.  The  spectrum  is  that 
for  a  single  point  on  the  ocean  surface  assuming  that  the  wind  has  been  blowing  with  an 
approximately  constant  speed  and  direction  for  a  period  sufficiently  long  for  the:  spec¬ 
trum  to  be  temporally  stable. 

This  spectrum  was  computed  from  data  assumed  to  be  representative  of  fully 
developed  seas,  generated  by  winds  of  20  to  40  knots.  Therefore,  the  data  on  wi  eh 
Eq.  (26)  is  based  are  not  representative  of  the  seas  analyzed  in  this  program— wh  ir 
3  to  21  knots  generating  seas  that  were  probably  only  partially  developed.  The  -  of 


r3?r 


■\\ 
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the  spectrum,  however,  which  has  not  been.shpwn  to  be  invalid  at  lower  wind  speeds, 
yields  results  consistent  with-,  visual  observations.  For  the  sea  conditions  encountered  in 
this  program,  an  experienced  observer  can.  estimate  the  wind  speed  from  whitecap  cover¬ 
age.  In  all  data  secs  .ftowing  wfeitecapping  the  wind  speed  estimated  from- the  appearance 
of  the  negatives  agreed  with  that  inferred  by  the  type  of  analysis  used  :here  to  better  than 
±  2  knots.  In  addition  to  this  type  of  analysis  being  more  accurate  at  the  lower  wind 
speeds,  the  coirAation  between  the  appearance  of  the  ocean  surface  and  the  inferred  wind 
speed  is  ve’  y  •  ood.  It  is  therefore  considered  that  the  use  of  the  Pferson-Moskowitz 
spectrum  m  .his  program  is  valid. 


t'y 

3 


"  '  following,  development  will  show  how  a  height  spectrum  may  be  transformed 
inV  ope  spectrum.  As  a  definition,  consider  any  quantity  of  the  form  b.  c )  to 
i>  spectrum  of  the  measurable  a  as  a  function  cf  b  and  c.  Since  both  a  one-dimen- 
d  height  spectrum  (as  in  Eq.  (26))  ana  a  two-dimensional  height  spectrum  (two- 
a.u'ensional  as;  in.  the  case  of  an  optical  spectrum)  of  the  same  surface  must  yield  the 
same  rn«s  wave  height  Oh, 


t2  = 


•* 

4>(h; 


u)  doi 


f  f2ff 

|  <b(h;  k,  a)k  dhda 

4)  ^0 


■  1 A»- 


a)k  dk  , 


where 


r Zn  , 

A(k,  a)  -  I  $(h;  k,  a)  da  (fe  =  ~ 


(27) 

(28) 

(29) 

(30) 


and 


a-*-**  . 


(31) 


Therefore 


<I>(ft;  to)  do) -=  A(k,  a)k  dk 


(32) 


For  the  wavelengths  being  analyzed  to2  =  gk.  Hence 

-2lf 


<l>[h;  io)  = 


4>(/i;  k,  a)  da 


(33) 


Since 


<&(h;k,a)  =  k-2']>{<P;k,a)  , 


(34) 


(43) 
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Knowing  the  form  that  I'fQi,  &w)  should  have,  the  wind  speed  V  may  be  determined 
by  fitting  the  two  curves,  using  either  Eq.  (41)  or  Eq.  (43).  The  following  discussion  will 
consider  only  the  Pierson-Moskowitz  form,  Eq.  (41).  Equation  (41)  may  be  rewritten  as 

*2 1}{k,^w)cce^s2/v4k2  (44) 


-."here  the  symbol  means  “varies  as.”  This  is  the  form  actually  used  in  fitting  the 
two  curves  to  determine  V.  The  exponential  rises  sharply  as  a  function  of  ft  and  then 
bends  to  a  slowly  increasing  function  for  large  ft.  The  bend,  or  knee,  of  the  curve  is  used 
to  fit  the  two  sides  of  Eq.  (44). 

The  data  used  by  Pierson  and  Moskowitz  to  obtain  the  wave-height  spectrum  were 
analyzed- separately  by  Moskowitz  [5]  to  yield  the  significant  wave  height  as  a  function 
of  wind  speed.  The  functional  relationships  between  the  significant  wave  height,  the  rms 
wave  height  o/,,  and  the  wind  speed  V  (Eq.  (1)  and  (10)  of 'Ref.  5)  are  such  that 

oh  =  0.00455  V2 

in  units  of  feet  if  V  is  given  in  knots. 


The  rms  slope  of  the  ocean  surface  is  also  a  function  of  wind  speed.  The  rms  slope 
in  the  upwind  direction,  the  rms  slope  in  the  crosswind  direction,  and  the  total  rms  slope 
have  been  determined  empirically  by  Cox  and  Munk  [6]  to  be 

rms  slope  (upwind)  =  0.1  J  0.163  V,  (46) 

rms  slope  (crosswind)  =  0.1  J 0.3  +  0.099  V,  (47) 


rms  slope  (total) 


0.3  +  0.264  V 


in  units  of  radians  if  V  is  given  in  knots.  The  wind  speeds  used  to  develop  the  proceed¬ 
ing  equations  were  measured  at  42  feet  above  the  surface.  The  difference  in  wind  speed 
however  between  64  feet  and  42  feet  above  the  surface  is  small.  In  addition,  the  rms 
slope  is  a  weak  function  of  wind  speed.  The  use  of  the  wind  speed  inferred  from  Eq.  (25) 
introduces  a  very  slight  error  that  is  within  the  experimental  error  of  this  type  of  analysis. 


The  Boeing  Company,  which  is  using  the  ground  truth  surface  descriptors  supplied 
by  NRL  to  correlate  electromagnetic  scattering  with  predictions  based  on  the  Kirchhoff 
approximation,  has  requested  that  rms  slopes  for  “slick  surfaces”  be  included  in  this  re¬ 
port.  The  Kirchhoff  approximation  assumes  the  radius  of  curvature  of  the  scattering 
facets,  and  hence  the  water  wavelengths  considered,  to  be  much  larger  than  the  electro¬ 
magnetic  wavelength.  The  electromagnetic  wavelength  at  L  band  is  about  8  inches;  thus 
it  may  be  necessary  to  consider  only  the  rms  slope  due  to  longer  wavelengths.  Rms  slope 

equations  resulting  from  data  taken  over  oil  slicks  have  also  been  given  by  Cox  and  Munk 

[6] .  The  oil  slicks  effectively  eliminated  wavelengths  shorter  than  1  foot,  yielding 

rms  slope  (upwind)  =0.1  J  0.5  +  0.040  ,  (49) 

rms  slope  (crosswind)  =  0.1  J  0.3  +  0.043  V ,  (50) 
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in  units  of  radians  if  V  is  given  in  knots.  The  values  reported  hereliave  been  calculated 
from  Eq.  (48)  and  (51). 

Recent  radar  data  and  a  few  spectral  forms  suggest  values  for  the  rms  slope  50% 
greater  than  those  given  by  Eq.  (48).  However  Eq.  (48)  and  those  indicating  higher  rms 
slopes  have  been  used  successfully  both  experimentally  and-in  various  theoretical  models. 
It  is  therefore  expected  that  the  true  rms  slope  for  a  clean  surface  would  be  in  the  range 
from  those  figures  given  in  this  report  to  values  50%  higher.  It  is  not  known  whether 
a  like  adjustment-in  Eq.  (49)  through  (51)  would1  be  valid. 


The  total  rms  slope  associated  with  a  given  band  of  water  wavelengths  Xi  to  X2  may 
also  be  calculated  from  the  equation 


cX2 

[rrns  slope  (Xj,  X2)] 2  =  9.27  X  10“4  V  j  x~  'e~xdx , 

Jxi 


(52) 


where  the  rms  slope  is  in  radians  if  V  is  given  in  knots  and  where 


and 


x2  =  3.60 


(53) 

(54) 


in  which  X  is  given  in  feet.  Equation  (49)  may  be  derived  by  multiplying  Eq.  (42)  by  k2 
and  integrating  from  cu-j  to  cj2  using  102  =  gk,  x  -  2g2IV2w2,  and  k  =  27 r/X.  This 
equation  is  not  valid  for  narrow  wavelength  bands  below  1  or  2  feet.  When  Xj  =  0  and 
X2  =  °°  Eq.  (52)  yields  a  total  rms  slope 

rms  slope  (total)  =  0.1  J 0.164  V  .  (55) 

This  value  is  low  compared  to  that  given  by  Cox  and  Munk,  which  itself  may  be  low. 

Any  value  calculated  using  Eq.  (52)  should  therefore  be  increased  by  an  appropriate, 
wind-speed  and  bandwidth  dependent,  factor.  This  factor  is  unknown  at  present.  Its 
value  may  be  appro.' imated  by  comparing  various  figures  given  by  Eq.  (52)  to  the  values 
given  by  Cox  and  Munk. 

The  wave-height  spectrum  depends  on  the  wind  speed  at  64  feet  and  may  be  given 
as  a  function  of  wave  number.  The  functional  relationship  is 

<l>(/i;  k)  =  0.0041  /r3e95/v4fc2 ,  (56) 

which  has  the  units  of  feet  cubed  when  k  is  given  in  inverse  feet  and  V  is  given  in  knots. 


COMPUTER  ANALYSIS 

There  are  three  circumstances  in  which  the  optical  diffraction  pattern  cannot  be 
properly  analyzed.  The  first  of  these  circumstances  occurs  when  the  knee  of  the  curve, 

Eq.  (44),  is  at  the  long  water  wavelengths  and  is, hidden  within  the  dc  term.  The  intensity 
distribution  beyond  the  knee  cannot  be  fitted  to  the  known  spectrum  with  any  degree  of 
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validity.  The  second  . circumstance  arises  when  the  diffraction  pattern  is  weak- and  most 
of  the  distribution  is  at  intensities  below  the  noise  level  ox  the  detector.  In  such  cases 
the  distribution  cannot  be  adequately  defined.  The  third  circumstance  occurs  when 
whitecaps  are  present.  They  appear  on  the  film  as  randomly  shaped  and  spaced  black 
patches  having  a  very  high  contrast  with  respect  to  the  waves.  The  whitecaps  therefore 
produce  a  random  diffraction  pattern  which,  due  to  the  higher  contrast,  is  sufficiently 
high  in  intensityto  partially  mask  that  due  to  the  surface  slope  angles.  The  diffraction 
pattern  due  to  the  whitecaps  cannot  be  included: in  the  analysis  because  the  form  of  the 
slope  spectrum  to  which  the  patterns-sre  fitted  was  derived  from  a  wave-height  spectrum. 
It  does  notinclude  the  random  and  wildly  fluctuating- slopes  present  in  a  wliitecap.  If 
one  attemptsto  optically  analyze  an  xsh  of  i  he  surface  shov.ing:no  whiteeapping,  one 
must  exclude  the  slightly  rougher  areas  of  the  ocean  surface  in  the  vicinity  of  the  white- 
caps  that  are  included  in  the-derived  slope  spectrum.  The  surface  descriptors  so  deter¬ 
mined  would  be  those  of  a  much  calmer  sea. 


Those  sets  of  photography  falling  into  one  or  more  of  the  above  categories  were 
analyzed  by  computer.  Each  set  was  inspected  visually  and/or  qualitatively  by  optical 
analysis  to  determine  the  direction:  of  the  wind.  A  microdensitometer  was  then  used  to 
obtain  a  trace  of  optical  density  versus  distance  bn  the  film  from  near  to.far  range.  A 
single  trace  in  the  direction  of  the  wind  was  obtained  for  each  negative.  The  aperture  at 
the  negative  was  a  slit  0.05  mm  wide  in  the.  direction  of  the  trace  and  0.60  mm  long. 

The  slit  was  long  enough  to  integrate  out  the  slopes  in  directions, other  than  the  direction 
of  the  wind  and  narrow  enough  to  reproduce  the  desired-slopes  with  adequate  fidelity. 
Individual  whitecaps  could  then  be  removed  by  drawing  dummy  slope  profiles  under  each 
whitecap.  The  profiles  on  each  side  of  the  whitecap  would  remain  to  contribute  to  the 
spectrum.  The  traces  were  then  digitized.  Each  trace  covered  approximately  2/3  the 
length  of  the  negative.  This  was.  done  to  adequately  sample  the  long  wavelengths. 


Since  the  length  of  the  trace,  approximately  25  mm,  was  more  than  double  the 
linear  dimension  of  the  area  used  in  optical  analysis,  nominally  11  mm,  three  steps  had 
to  be  taken  to  reduce  distortions  in  the  subsequently  computed  spectrum.  The  overall 
slope  of  the  trace,  D0(x,  y)  in  Eq.  (10),  and  the  average  value  were  removed.  These  two 
steps  greatly  reduce  the  spectral  “noise”  at  low  frequencies  and  produce  the  equivalent  of 
a  much  tighter  dc  term  in  optical  analysis.  Perspective,  which  shortens  wavelengths  at  far 
range  compared  to  those  at  near  range,  was  also  removed.  This  is  necessary  when  analyz¬ 
ing  a  long  trace  to  prevent  overintegration  in  the  spectrum. 


The  density  differences  embodied  in  the  final  traces  were  small;  consequently  the 
traces  were  equivalent  to  those  of  <P(x,  y)  versus  distance  on  the  mean  ocean  surface. 

Each  trace  was  then  analyzed  by  means  of  a  Fast  Fourier  Transform  program  to  yield  a 
curve  equivalent  to  the  intensity  distribution  arrived  at  in  optical  analysis,  Eq.  (24).  Since 
the  individual  transforms  were  representative  of  only  single  traces  through  a  given  negative, 
the  curves  fluctuated  wildly  as  a  function  of  fe(=  2?r/X)  and  could  not  be  fitted  reliably  to 
a  known  spectrum.  Consequently  all  of  the  spectra  obtained  for  a  given  set  of  photo¬ 
graphs  had  to  be  averaged  to  yield  a  single  spectrum.  This  spectrum  was  then  fitted  to 
the  known  spectrum  to  determine  the  required  surface  descriptors. 


MINOR  SPECTRAL  DISTORTIONS 


Several  minor  distortions  in  the  experimental  spectrum  and  ways  to  avoid  their 
effects  will  now  be  discussed.  Most  of  these  arise  from  the  manner  in  which  the  sea  is 
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on°fh?^?hedV  bussed  from, the  viewpoint  of  optical  analysis  -the  effect 

on  the  computer  analysis  being  equivalent.  y  ’  ne  enect 

The  overall  density  gradient  on  the  negative  Dn(x  v*  in  Fn  n  m 
of  the  sensitivitv  fun pf inn  tj  •  n  *  O'  *  the  variation 

A  •  function,  U0(xt  y )  in  Eq.  (11),  have  been -discussed  in  the  section  on 

Optical  Analysis  and  have  been  shown  to  place  an.  inconsequential  limit  on  one’s  abilifcv 

,  ^fc’  !  !l  h  3150  ^cou-^iuent,  U  in  practice,  on  the  area  of  the  detector 

Neitte  °f  thKe  tw0  f“"clions  *-» 

dist^1Ce  f5T  thC  °rigin  in  thc  °Ptical  diffraction  pattern  is  proportional  to  the 

sr  "5  tecorded  °,n  the  fx  z\£ 

Sth  luTth^r  “nr“PS”dt°  Wgher frequency,  shcu.rt,  waves.  The 
he  :a"  SKVS*  ls  recorded  on  the  negative  thus  introduces  a 
distortion  into  the  spectrum.  It  the  same  wavelength,  were  present  everywhere  on  the 

?Th“Ve,  the  lK,VeS  at  •«' Saga  -oSKSS *gette^id 

ranS  ^£52*  ^hJ,fart.!?*  hT  the  0rigin  than  'v°uld  the  530)6  waves  at  nearr 
is  “smeared»1n  the  direcS  ”  °  trantsfo5m  pIane  corresponding  to  a  given  wave.3ngth 

ir0m  ne“  to  »  ra°gc-  ^  therefore  introduces  a  slight 
mtcgraumg  iaCtor  into  me  spectrum.  However,  it  is- usually  necessxirv  to  in fawrft-  fh- 
spectrum an  even  greater  degree,  through  the  use  of  a 

a  reasonably  smooth  curve  to  fit  to  the  known  spectral  form  *  ’  0Wa” 

rave^ri^is^sU^tl^cfetortei^to  lavw^slopes^i^tedtcwM^^^ramCTd^^ThJdistortion 

Sn^tSr^n8ity  3dded  *  ^  «**■?«“"  Po^on  of  the  spectrum 

■ii^-«Steeper  than  ROrmal  siopos  that  are  associated  with  the  shorter  wavelengths 
he  d  fr?Ct  “  °T  °f  light  energy  into  the  highT^Tporta  of 

Z  Pa“e?'-  ™*  °£  courso  on  "better  the  short  wavefengThs  ie 

f  °n  jj9  negative.  If  these  short  waves  are  present  but  not  resolved  thev  mav 

Sglcs'  A °fa«t  of  thfsurfa^  WgWf ’"'t11*  !**«<»>  by  distorting  the  recorded  slope 
angles,  a  tacet  of  the  surface  covered  by  unresolved  waves  is  recorded  at  a  densitv 

nd.cat,ve  of  a  facet  sloped  slightly  more  to  ward  the  car-era.  ThTsSp  dopes tX  in 

frequently  occur  on  the  longer  wavelengths,  and  are  not  adequately  represented  bv  the 

MornpatteTnyFllm  T? Z  ta‘e““?  ”  the  bigh-ucquency  portion  of  the  dif 

naction  pattern.  Film  gia^  too  increases  the  intensity  in  the  high-frequenev  oortinn  of 

trtm^l0nr^  F,lm  grain  has  essentially  the  same  effect  as  a^arge  mLber  of 
a  ly  uniformly  dispersed,  extremely  small  whitecaps.  The  diffraction  pattern  of  the 

,  1S  randfom>  since  the  grains  are  of  necessity  smaller  than  the  finest  detail 

e  anaiyzed,  it  contributes  to  the  spectrum  only  at  very  high  spatial  frequencies  The 

ZT6  Ifw the  atmospnere  when  the- photographs,  are  'aken  has  the  opposite 
effect  on  the  high-frequency  portion  of  the  diffraction  pattern.  It  reduces  ^resolution 
on  e  negative  and  smooths  the  a  .gularity  of  the  randomly  deformed  water  surface  Th« 
intensity  in  the  high-frequency  P<  rtion  is  thus  reduced  by  the  presence  of  haze 

in  the  preceeding  Paragraph  make  the  high-frequency 
portion  of  the  experimental  spectrum  somewhat  ambiguous.  This  (an  be  avoided  bv 

placing  the  emphasis  of  the  analysis  on  the  low-frequency  portion.  HuS  thetee  of  the 
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curve,  Eq.  (44),  falls  within  the  low-frequency  portion  of  the  spectrum  in  the  vast  majority 
of  cases,  the  ambiguity  is  of  little  consequence  in  this  type  of  analysis; 

DATA  ANALYSIS 

Ah  inventory  of  test  film  and  photographs  taken  by  FAA  and  submitted  for  analysis  is 
given  in  Table  1.  The  first  nine  rolls  of  film  were  obtained  in  testing  the  camera  system  to  be 
used  and  in  gaining  experience  in  the  techniques  of  photography,  necessary  for  Sea  Photo 
Analysis.  Twenty  sets  of  data  film  and  three  rolls  of  test  film  were  submitted  after  the  pre¬ 
liminary  test  period.  Four  of  the  data  sets  were  not  suitable  for  analysis,  as  explained  below. 
A  representative  aerial  photograph  from  each  of  the  sixteen  data  sets  analyzed  is  shown  in 
Figs.  7a  through  7p,  The  photographs  were  taken  from  an  altitude:  of  1000  feet,  showing  an 
area  roughly  600  feet  by  1600  feet  at  a  mean  depression  angle  of  35°  (40°  in  Figs.  7m 
through,  7p). 


Table  1 

Inventory-of  Data  Film 


Comments* 


4-16-71  6 

7  - 

4- 30-71 

5- 1-71  1 

2 

5-24-71  1, : 

5-25-71  2, : 

5-26-71  4, ! 

7-12-71 
7-26-71 

9-1-71  1 

9-2-71  2 

9-3-71  3 

9- 22-71 

10- 13-71  1 

10-14-71  2 

10-15-71  3 

3-4-72  1 

'  3-7-72  2 

3- 8-72  3 

4- 20-72  1, ! 

4-21-72  3,. 

4-22-72  5,  i 

4-24-72  7, : 

4- 24-72  9 

5- 17-72  1 

5-18-72 _  2 

*The  comments  are  explained  in  the  test. 


Resolution  test 
Overexposed;  smrglitter 
Overexposed,  sun  glitter 
Overexposed 

Overexposed,  camera  tests 
Camera  tests 
Overexposed,  sun  glitter 
Camera  tests 
Test  film 

Group  4  analyzed  optically 
Group  7  analyzed  by  computer 
Group  4  analyzed  by  computer 
Group  5  analyzed  by  computer 
Camera  tests 
Camera  tests 
Degraded  imagery 
Group  3  analyzed  optically 
Degraded  imagery 
Group  6  analyzed  optically 
Group,3  analyzed  optically 
Group  3  analyzed  optically 
Group  3  analyzed  optically 
Group  6  analyzed  by  computer 
Group  5  analyzed  by  computer 
Group  6  analyzed  optically 
Group  2  analyzed  by  computer 
Group  2  analyzed  optically 
Group  1-  analyzed  by  computer 
Group  1  analyzed  optically 
Camera  calibration 
Wrong  polarization 
Wrong  polarization 


NRL  REPORT7574 


22  ROGER  O.  PILON 

As  rioted  in  Table  l,the-photographsin  three  of  the  preliminary  data  sets  showed 
sun  glitter  patterns.  Sun  glitter  is  the  specular  reflection  of  the  sun  directly  toward  the 
camera.  On  the  negative  it  appears  as  black  spots  distributed  randomly  around  a  central  (  ; 

point.  Ingeneralthe  spots  decrease  in  size  and  are  spaced  farther  apart  with, increasing 
"distance  from -the  central,  point.  Since  the  distribution  is  random  from  frame  to  frame,  ; 

the  diffraction  patterns  are  likewise  random.  Also,,  since  the  contrast  between  tha  very  ! 

black  spots  and.  the  much  lighter  surroundnig  image  of  the  waves  is  much  greater  than  ' 

that  due  to  wave  slopes,  the  intensity  of  the  extraneous  diffraction  pattern  due-to  the  • 

glitter  partially  masks  that  due  to  the  waves.  These  two  factors  would  make  it  very  ; 

difficult  to  extract  the  desired  spectrum.  In- addition,  those  slonos  reflecting  sunlight  j 

directly  toward*  the  camera  axe  effectively  removed  from  the  slope  analysis,  thus  distort-  ! 

ing  the  spectrum  somewhat.  For  these  reasons  it  is  necessary  to  obtain  photographs  of  | 

the  ocean  surface  that  do  not  show  glitter  patterns.  This  still  leaves  approximately  270°  » 

of  the  360°  of  azimuthal  directions  available  for  photography.. 

The  photographs  obtained  on  1  and  3  September  1971  were  degraded  by  ccndensa-  j 

tion  and  hydraulic  oil  on  the  window  through  which  the  ocean  surface  was  photographed.  j 

Foreign  fluids  on  the  window  effectively  smear  the  images  of  the  shorter  wavelengths,  re-  j 

moving  them  from  the  analysis,  and  coiild  introduce  density  variations  on  the  film  that  j 

would  distort  the  long  wavelength  analysis.  Thus  such  data  sets  are  not  suitable  for  i 

analysis. 

The  photographs  of  17  and  18  May  1972  were  obtained  with  the  polarizer  *>Iaeed 
over  the  camera  lens  in  the  wrong  orientation.  The  surface  was -therefore  photographed 
using  vertically  polarized  light  instead  of  horizontally  polarized  light.  The  intensity  re¬ 
flection  coefficient  curve  corresponding  to  Fig.  3  for  horizontally  polarized  light  would 
have  the  same  values  at  0°  and  90°,  but  would  go  to  zero  at  Brewster’s  angle,  53°. 

Therefore  R(f J),  maintaining  a  mean  depression  angle  of  40°,  would  be  double  valued.  In 
addition  the  mean  light  level  and  both  the  magnitude  and  directionality  of  the  wave-slope 
sensitivity  would  change  as  a  function  of  position  on  the  film.  The  theory  required  to 
analyze  such  data  sets  has  not  been  developed. 

Eight  of  the  data  sets  were  analyzed  optically,  as  noted  in  Table  1.  The  photo¬ 
graphs  obtained  on  1  May,  22  September,  and  13,  14,  and  15  October  19.71  contained 
an  insufficient  number  of  analyzable  frames  to  warrant  temporal  plots  of  the  ocean 
descriptors.  The  minimum,  average  and  maximum  values  of  the  descriptors  obtained  from  ‘ 

the  few  usable,  frames  are  given. in  Table  2.  The  five  sets  of  descriptors  inferred  from  the  j 

1  May  1971  photography  is  given  in  Table  3. 

! 

The  surface  descriptors  inferred  from  the  three,  remaining  data  sets  that  could  be  j 

analyzed  optically,  8  March  1972  and  21  and  24  April  1972,  are  given  in  Table  2  and  i 

are  plotted  as  a  function  of  time  and  distance  in  Figs.  8  through  16.  Error  bars  are  J 

shown  in  Figs.  8  and  9.  They  represent  the  expected  error  associated  with  only  one -step  ( 

in  theanalysis  procedure,  namely  that  of  filing  the  intensity  distribution  of  the  diffrac-  } 

tion  pattern  to  the  accepted  spectral  form,  and  are  indicative  of  the  corresponding  j 

expected  error  in  the  other  sets  reported  here.  They  are  shown  as  an  indication  of  the  J 

agreement  between  the  ocean  surface  to  be  described  and  the  spectral  form  used  to  > 

obtain  rms  wave  height  and  rms  slope. 
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Table2 

Inferred  Ocean  Surface  Descriptors 


Averaged 
Heading 
of  Waves 
(deg) 


Wind  Speed 
at  64  ft* 
(knots) 


RMS 

Wave  Height* 
(ft) 


RMS  Slope*  (radians) 


Clean 

Surface 


5-1-71 


5t24t71 

5-25-71 

5-26-71 

9-2-71 

9-22-71 


10-13-71  1821  180 


10-14-71  1827  135 


10-15-71 


3-4-72 

3-7-72 

3-8-72 


4-20-72 

4-21-72 


4-22-72 

4-24-72 


0.050 

0.069 

0.084 

2.1 

0.19 

0.70 

0.98 

0.19 

0.22 

0.26 

0.29 

0.52 

1.0 

0.39 

0.46 

0.57 

0.32 

0.46 

0.61 

1.5 

0.96 

0.17 

0.22 

0.33 

0.93 

0.20 

0.22 

0.23 

1.0 

0.12 

0.15 

0.17 


♦When  three  values  are  given  for  any  one  quantity  they  represent  the  minimum,  average,  and  maximum 
values  for  the  data  set. 

fRms  slope  due  to  wavelengths  longer  than  1  foot. 


Computer  analysis  was  required  on  -the  eight  remaining  data  sets,,  24,  25,  and  26 
May  1971,  2  September  1971,  4,  and  7  March  1972,  and  20  and  22.  April  1972;  Con¬ 
sequently,  only  the  average  values  for  wind  speed,  rms  wave,  height,  and  rms  slope  (clean 
and  slick  surface)  are  given  in  Table  2. 
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Table  3 

Comparison  of  Descriptors  InferredFror.i  the I'ierson-Moskow.itz  Spectrum 
and  From  the  Neumann  Spectrum* 


Time 

(GMT); 

Windspeed  at  64  ft 
(knots) 

-  — “I 

RMS  Wave  Height 
(ft) 

RMS  Slope  (radians)  j 

Clean  Surface 

Slick  Surface 

Pierson- 

Moskowitz 

Neumann 

Pierson- 

Moskowitz 

Neumann 

Pierson- 

Moskowitz 

Neumann 

Pierson- 

Mcskowit/. 

' 

Neumann 

18:11:59 

4.0 

4.3 

0.07$ 

0.084 

0.12 

0.12 

0.11 

0.11 

18:12:01 

3.8 

4.0 

0.066 

0.073 

0.11 

0.12 

0.11 

0.11 

18:12:03 

4.3 

4;5 

6.084 

0.092  - 

0.12 

0.12 

0.11 

0.11 

18:12:05 

4.0 

4.3 

0.073 

0.084 

0.12 

0.12 

0.11 

0.11 

18:12:07 

4.0 

4.3 

o:d73 

0.084 

0.12 

0.12 

C.ll 

0;11 

16:13:04 

3.3 

3.3 

0.050 

0.050 

0.11 

0.11 

0.10 

0.10 

♦Obtained  from  the  1  May  1971  Photographs  (Roll  2,  Group  4).  The  Averaged  Heading  of  the  Waves 
was  20°. 
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Fig.  8— The  inferred  wind  speed  as  a  fUftp-fan  of  time  and  distance  for  8  March  1972 


The  wind  speeds  inferred  from  Eq.  (44)  appear  to  be  quite  accurate,  based  on  the  fact 
that  an  experienced-observer  can  estimate  the  wind  speed  quite  well  when  whitecaps  are 
present.  In  all  data  sets  showing  whitecapping  the  estimate  agreed  with  that  derived  from 
Eq.  (44)  to  better  than  ±2  knots.  In  addition  to  this  type  of  analysis  being  more  accurate  at 
the  lower  windspeeds,  the  correlation  between  the  appearance  of  the  ocean  surface  and  the 
inferred  wind  speed  is  very  good.  The  values  computed  for  the  rms  wave  height  and  rms 


(radiant) 
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0  1  2  3  4  5  6 

DISTANCE  (nautical  miles) 

Fig.  15— The  inferred  rms  wave  height  as  a  function  of  time  and  distance  for  24  April  1972 


0  1  2  3  4  5  6 


DISTANCE  (nautical  miles) 

Fig.  1G— The  inferred  rms  slope  as  a  function  of  time  and  distance  for  24  April  1972 

Pertinent  facts  concerning  the  photographs  have  been  explained  and  the  type  or 
analysis  used  for  each  data  set  has  been  listed.  The  time  of  day,  averaged  heading  of 
waves,  wind  speed,  rms  wave  height  and  rms  slope  for  clean  and  slick  surface,  have  been 
determined  for  each  analyzable  data  set,  with  minimum,  average,  and  maximum  values 
given  for  those  sets  that  could  be  analyzed  optically.  Root-mean-square  slopes  due  to 
different  water  wavelength  bands  and  wave-height  spectra  may  be  calculated  using  the 
equations  given  in  this  report.  The  parameters  inferred  from  the  analysis  appear  to 
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describe  very  well  the  ocean  surface  at  the  time  and  place  at  which  the  photography  was 
taken. 
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APPENDIX 


HINDCAST  DATA  AND  COMPARISON  WITH  INFERRED  DESCRIPTORS 

The  techniques  of  Sea  Photo  Analysis  have  been; used  to  obtain  the  required  descrip¬ 
tors  of  the  wind-generated  seas  encountered  in  this  program.  The  roughness  known  as 
sea,  however,  is  but  one  of  the  two  components  of  the  ocean  surface  affecting  the  scatter¬ 
ing  of  electromagnetic  energy.  The  other  component  is  swell.  Therefore  the  Naval  Re¬ 
search  Laboratory  recommended  that  the  Federal  Aviation  Administration  obtain  hindcast 
data  on  swell.  Such  data  were  obtained  from  the  Naval  Oceanographic  Office  and  are 
summarized  in  Table  Al.  The  data  give  the  direction  from  which  the  swell  was  propagat¬ 
ing,  the  wave  height  measured  from  crest  to  trough,  and  the  period  of  the  swell.  The 
wavelength  L  of  the  swell  (assumed  to  be- sinusoidal)  may  be  estimated  from  [Al] 

L  =  5.12  T2  ,  (Al) 

where  L  has  the  units  of  feet  when  the  period  T  is  given  in  seconds. 

Table  Al  also  includes  hindcast  data  for  the  sea 'by  the  Naval  Oceanogtaphic  Office 
(NOO)  and  for  a  combined  sea  and  swell  by  the  Canadian  Forces  Meteorological  and 
Oceanographic  Center,  Atlantic  (CFMOCA).  Three  differences  exist  in  the  manner  of 
reporting  the  sea  (or  combined  sea  and  swell)  hindcasts  and  the  inferred  descriptors  given 
in  this  report.  First,  the  time  and  location  for  which  the  hindcast  is  given  is  not  always 
the  same  as  the  time  and  location  at  which  the  aerial  photographs  were  taken.  The 
positions  for  the  hindcasts  differ  by  as  much  as  250  nautical  miles  from  those  at  which 
the  photographs  were  taken.  Depending  on  frontal  positions  and  wind  field  perturbations 
in  the  area,  this  could  account  for  some  of  the  discrepancy  between  the  hindcast  and  the 
inferred  data.  Second,  the  hindcasts  give  the  direction  from  which  the  wind  was  blowing, 
generating  (sea)  waves  heading  in  the  opposite  direction.  Therefore,  the  directions  given 
in  the  sea  (or  sea  and  swell)  hindcasts  differ  by  180  degrees  from  those  given  in  the  main 
text  (tables  2  and  3).  Third,  the  hindcasts  list  a  significant  wave  height,  H ,  ,3  (the 
average  of  the  highest  1/3  of  all  wave  heights  measured  crest  to  trough),  which  is  related 
to  the  rms  wave  height  given  in  this  report  by  [A2] 

#1/3=40,,.  (A2) 

The  inferred  descriptors  given  in  this  report  are  compared  to  the  sea  (sea  and  swell) 
hindcasts  in  Table  A2.  A  linear  interpolation  in  time  was  used  to  adjust  the  NOO  hind¬ 
casts  to  the  time  the  photographs  were  taken.  The  differences  in  wave  headings  were 
determined  by  adding  180°  to  the  directions  given  in  the  hindcasts  and  then  subtracting 
that  value  from  the  Averaged  Heading  of  Waves  given  in  Table  2.  The  headings  given  in 
Table  2  were  determined  from  the  photographs  and  spectra  to  better  than  ±3°.  There¬ 
fore,  the  differences  in  wave  headings  given  in  Table  A2  are  mostly  errors  in  the  hindcast 
data  (except  in  those  cases  where  the  time  and/or  location  for  the  hindcast  is  sufficiently 
removed  from  that  at  which  the  photographs  were  taken  to  be  in.  a  different  wind  field). 
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Table -A2 

Comparison  of  Hindcast  Data-With  Inferred  D  escrip  tors 


Date 

- 

Difference  in  Wave. 

Heading  (deg) 

RMS  'Waye-HeightRatio  | 

NRL-NQO 

NRL- CFMOCA 

NRL/NOO 

NRL/CFMOCA 

5-1-71 

-25 

0:055 

'/ —  ; 

5-24-71 

+16 

— . 

6:51 

- 

5-25-71 

+4i 

— 

6.50 

■ 

5-26-71 

-25 

— 

6:67 

rr  ' 

9-2-71 

tli 

— 

0.33 

— 

9-22-71 

-88 

1.8 

046: 

10-13-71 

-45 

-90 

0:82 

0,51 

1044^71 

+35 

— 

0.43 

1045-71 

+90 

-60 

0'7  T 

0:56 

3-4-72 

+96 

- 

2,2 

— 

3-7-72 

+17 

+47 

1.0 

0.96 

3-8-72 

+68 

-60 

0:40 

0:95 

4-20-72 

+33 

-07 

0.26 

2:2 

4-21-72 

-35 

+70 

.6:19 

0.078 

4-22-72 

+19 

+03 

0.54 

0:92 

4-24-72 

-06 

+80. 

0:17 

0.11. 

-»  *  •  - 

The  rms  wave-height  ratio  is  the  inferred  rms  wave  height  divided  by  the  hindcasted 
rms  wave  height  for  the  sea  (W  the  significant  wave  height  of  the  sea).  The  wave  heights 
hindcasted  by  the  Naval  Oceanographic  Office  for.  swell  were  subtracted  from  the  hind- 
casts  by  the  Canadian  Forces  Meteorological  and  Oceanographic  Center,  Atlantic  to  yield 
an  effective  significant  wave  height  for  the  sea.  These  effective  values  were  then  used  to 
calculate  the  rms  wave-height  ratios  for  the  NRL/CFMOCA  data.  In  seven  of  the  com¬ 
parisons,  1  May  1971  (NOO),  22  September  1971  (CFMOCA),  44  October  1971 
(CFMOCA),  and  21  and  24  April  1972  (NOO  and- CFMOCA),  visual  inspection  of  the 
photographs  reveals  that  the  hindcasts  do  not  describe  the  photographed  conditions.  The 
effective  wind  speeds  required  to  generate  the  significant  wave  heights  hindcasted  for  these 
sets  may  be  calculated  using  Eqs.  (A2)  and  (45).  The  wind  speeds  so  calculated  range 
from  14  to  25  knots.  Since  whitecapping  begins  at  wind  speeds  of  about  10  to  12  knot-s, 
the  photographs  for  the  above  sets  should  show  whitecapping.  They  do  not.  (The  imper¬ 
fections  in  the  unretouched  examples  in  Fig.  7,  such  as  film  scratches  and  dust  spots, 
should  not  be  confused  with  whitecapping  as  shown  in  Fig.  7b,  7d  (very  few),  7e,  7g,  7j, 
7k,  7m,  and  7o.)  The  effective  wind  speeds  for -these  sets  are  therefore  less  than  12 
knots,  as  indicated  by  inferred  wind  speeds  in  Table  2. 

The  effective  wind  speeds  required  to  generate  the  significant  wave  heights  hindcasted 
for  2  September  1971  and  20  April  1972  (NOO)  are  26  and  28  knots  respectively.  Therer 
fore,  Figs,  7e  and  7m  should  show  more  whitecapping  and  spray  than  Fig.  7b  (21.4  knots). 
Since  these  figui  3  show  significantly  less  whitecapping  and  no  spray,  the  hindcasts  are  not 
representative  of  these  data  sets. 

In  the  remaining  comparisons  the  rms  wave-height  ratio  varies  between  0.40  and  2.2. 
In  view  of  the  foregoing,  such  differences  in  data  taken  on  site  and  that  hindcasted  from 
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macroscopic  pressure  systems  and  wind  fields  over  the  open  ocean  should  be  considered 
very  good. 
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